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Abstract 

We prove solvability theorems for relaxed one-sided Lipschitz mul¬ 
tivalued mappings in Hilbert spaces and for composed mappings in 
the Gelfand triple setting. From these theorems, we deduce proper¬ 
ties of the inverses of such mappings and convergence properties of a 
numerical scheme for the solution of algebraic inclusions. 
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1 Introduction 

The relaxed one-sided Lipschitz property (see Definition [4] below) was first 
considered in [9], where it was identified as an important stability criterion 
for time-dependent differential inclusions. The behavior of general multival¬ 
ued mappings with negative relaxed one-sided Lipschitz constants was later 
studied in P33HH3- In particular, surjectivity of the mappings and therefore 
solvability of the corresponding algebraic inclusions was shown by consider¬ 
ing the flow of the differential inclusions from [9]. However, no information 
on the localization of the solutions was given in these papers. 

For relaxed one-sided Lipschitz mappings in finite-dimensional spaces, 
the solvability theorem given in [3j specifies a ball in which a solution of the 
inclusion is contained. The radius of this ball depends on the norm of the 
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residual of the inclusion at the center point. This theorem guarantees that the 
implicit Euler scheme for stiff ordinary differential inclusions is well-defined 
and convergent on the infinite time interval, and it has recently been applied 
in [T3] to obtain a numerical method for the solution of the generalized Bolza 
problem. A refined solvability result presented in [5] and restated as Theorem 
[TUI below immediately gives rise to a numerical algorithm for the solution of 
algebraic inclusions. 

These solvability theorems are relevant for the following reason. For non¬ 
scalar mappings, it is currently unclear whether continuous and relaxed one¬ 
sided Lipschitz multivalued mappings possess parameterizations that are con¬ 
tinuous and one-sided Lipschitz with the same one-sided Lipschitz constant. 
Moreover, simple examples show that selections generated by metric projec¬ 
tion such as the minimal selection are not one-sided Lipschitz with the same 
constant as the multimap. It is therefore impossible to obtain precise solvabil¬ 
ity results by applying standard tools like topological fixed point theorems 
to selections or parameterizations of one-sided Lipschitz multifunctions. 

In the present paper we generalize the finite-dimensional solvability re¬ 
sult from [5] to infinite-dimensional Hilbert spaces, and we discuss implica¬ 
tions both in an abstract framework as well as in the context of a special 
class of systems of elliptic differential inclusions. After collecting definitions 
and preliminary tools in Section [2], we prove an abstract solvability result 
in an infinite-dimensional Hilbert space in Section [3] via an approach based 
on Galerkin approximations. The approach avoids strong compactness as¬ 
sumptions which are not satisfied in many applications. In Section [4l we 
reformulate the main result in the context of Gelfand triples and composed 
multivalued operators. In this setting, special care was taken to obtain op¬ 
timal estimates by considering a mixed scalar product adapted to the prop¬ 
erties of the individual operators. As a byproduct, the main result reveals 
certain aspects of the behavior of the inverses of relaxed one-sided Lipschitz 
mappings as detailed in Section [5] As in the finite-dimensional context, the 
solvability theorem gives rise to a numerical algorithm for the solution of re¬ 
laxed one-sided Lipschitz algebraic inclusions, which is analyzed in Section [6j 
In Section [TJ we discuss a system of elliptic differential inclusions where the 
assumptions of the Gelfand triple version of our main result are verified for 
suitable right-hand sides. Moreover, we test the numerical algorithm from 
Section in the context of this system. 
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2 Preliminaries 


In this section, we collect the necessary definitions and some elementary facts. 
Let ( X , || • ||x) be any real norrned vector space, and let (-, •) : X* xI->R 
denote the dual pairing. 

Definition 1. For x G X and nonempty subsets M,M' C X, we set 
dist x(x,M'):— inf ||x — x'\\x, 

x'eM' 

e x {M, M') := sup distx(x, M'), 

xeM 

||M|| :=e x (M,{0 X }), 

Proj x (x, M) := {x' G M : \\x — x'\\x = distx (x, M)}, 

B X {M,R ) ■= {y E X : dist x (y,M) < Ft}, 

B x (x, R) := {y G V : \\y - x\\ x < R}- 

The nonempty closed, bounded and convex subsets of X are denoted CBC(X), 
and the nonempty convex and compact subsets of X are denoted CC(X). 

Definition 2. a) The support function o* x : X xCBC(X*) —> R is defined 
by 

crx(x, A) := sup(<p, x) Vx G X, A C X*. 

y&A 

b) If X is a real Hilbert space with scalar product (• j -)x, then we define 
(7 x : X x CBC(X) -> R by 

Ux(x, A) := sup(|/, x) Vx G X, A C X. 

y eA 

Definition 3. Let (M, d ) be a metric space and Y a further norrned vector 
space. 

a) A set-valued mapping F : M — CBC{Y*) is called upper hemicontinu- 
ous (uhc) at x G M if for any sequence (xfc)fc 6 N C M with x& —> x we 
have 

lim sup (Tyiv, F(xk)) < erf (y, F(x)) for all v G Y. (1) 

k —^oo 

It is called uhc if it is uhc at any x G M . 

If M is weakly sequentially closed, then F is called compactly upper 
hemicontinuous (c-uhc) if condition (pQ) holds for any sequence (xfc)fc €]N C 
M with Xk —^ x G M . 
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b) If Y is a Hilbert space, then a set-valued mapping F : M —$■ CBC(Y ) 
is called upper hemicontinuous (uhc) at x G M if for any sequence 
(^fc)fcGN C M with Xk —> x we have 

limsup a x {v, F(xk)) < <Jx(v, F(x)) for all v G Y. (2) 

k—> oo 

It is called uhc if it is uhc at any x G M. 

If M is weakly sequentially closed, then F is called compactly upper 
hemicontinuous (c-uhc) if condition (J2]) holds for any sequence (xk)ke¥i C 
M with Xk — 1 x G M. 

c) A set-valued mapping F : M —)■ CC{Y ) is called upper semicontinuous 
(use) at x G M if for any sequence ( Xk)k£¥i C M with Xk —> x G M we 
have 

ev{Fx(xk), Fn(x)) —>• 0 as k ^ oo. (3) 

It is called use if it is use at any x G M. 

The following one-sided property is the central object of investigation in 
the present paper. 

Definition 4. A mapping F : X =4 is called l-relaxed one-sided Lipschitz 
with constant l G R ("or l-ROSL) if for any x,x' G X and y G F(x) i/iere 
exists y' G F(x') suc/i that 

(y' - y,x' - x) < l\\x' - x\\ 2 x . 

In Theorems [9l [TO] and TU the relaxed one-sided Lipschitz property will 
only be required relative to one point in the graph of F. In Sections [5] and 
[6] however, we will deal with mappings that are relaxed one-sided Lipschitz 
in the sense of Definition [4] 

Remark 5. The definition of an ROSL mapping with constant l > 0 is for¬ 
mally similar to that of a monotone mapping. Nevertheless, ROSL and mono¬ 
tone mappings have fundamentally different properties. Monotone mappings 
on Hilbert spaces are single-valued outside a set of first Baire category (see 
m, and the operator IFolT is onto and possesses a single-valued inverse for 
any monotone T and any a > 0 , which is the theoretical basis for the proximal 
point algorithm (see f73| /). In contrast, the 1 -ROSL mapping F : R —» CC(R) 
given by F(x) = x + [—1,+1] is set-valued on the whole space, and for all 
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a > 0 the inverse ofl + aF, given by (/ + ctF){x) = (1 — a)x+ [— a, a], is set 
valued as well. Similarly, ROSL mappings with constants l < 0, which are 
mainly considered in the present paper, look formally similar to mappings T 
with —T monotone but have fundamentally different properties. 

We also recall the following facts which are well known and easy to see. 

Lemma 6. Let X be a Hilbert space, and let x E X. 

a) If M C X is closed and convex, then Proj(x,M) is a single point. 

b) If M C X is weakly sequentially closed, then Proj (x,M) is nonempty. 

c) If M C X is closed, and if (x n ) n C X and x G X satisfy \\x n — x IU ->• o 
and distx(.'E n , M) —> 0 as n —>■ oo, then x E M. 

The following standard observations will also be used later on. 

Lemma 7. Let Y be a reflexive Banach space, Z a normed vector space 
and T G C(Y,Z). Then for every A G CBC(Y) we have T(A ) C CBC(Z). 
Consequently, every map F : M —>■ CBC(Y ) defined on an arbitrary set M 
gives rise to a map T o F : M —> CBC(Z) 

Proof. Let A G CBC(Y). Since T is linear and continuous, T(A) is convex 
and bounded. To see that T(A) is closed, we consider a sequence ( Zk)k in 
T(A) such that Zk —> * G Z as k —> oo. Choosing G A such that T[yf) = Zk 
for k E N, we obtain a bounded sequence (yk)k hi A. Since Y is reflexive, 
we may pass to a subsequence such that yk —^ y G Y as N' 3 k —> oo, and 
y G A by Mazur’s Theorem. Moreover, 

T{y) = w-lini/^oo Ty k = w-lim^oo z k = z 

and thus z G T(A). Hence T(A) is closed. □ 

Lemma 8. Let X be a reflexive Banach space, and let M',M" G CBC(X). 
Then 

M' + M" := {m' + m" : m E M , m" E M"} E CBC(X). 

Proof. It is easy to check that M' + M" is bounded and convex. We show 
that M' + M" is closed. Let (m n ) n C M' + M" be any sequence with 
lim rwoo m n = m E X. Then there exist ( m' n ) n C M' and (m")„ C M" such 
that m n = m' n +m'f for all n. By the Banach-Alaoglu theorem, we have m' n —^ 
ml E M' along a subsequence, so that m" = m n — ml n —^ m — ml. By Mazur’s 
lemma, m" := m — ml E M". Therefore, m = m' + m" E M' + M". □ 
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3 An infinite-dimensional solvability theorem 

Let V be a separable Hilbert space with scalar product (-, -)y, associated 
norm || • ||y. The main result of this section is the following solvability 
theorem. 

Theorem 9. Let x, y G V, R > 0 and l < 0, and let F : B v (x, R ) —> CBC(V) 
be a multivalued mapping. 

a) Let F be bounded and c-uhc. If there exists some y G F(x) such that 
||y — y\\v < — IR and 

Vx G B v (x, R) 3 y G F[x) : (y — y,x — x)y < l\\x — x\\y, 
then there exists some 

x€B v (x c ,-±\\y-y\\ v ) with x c = x + ±(y - y), 
satisfying y G Fix). 

b) If F admits a modulus of continuity relative to x in the sense that 

ev(F(x), F(x)) < oa{\\x — x\\v) 

for all x G B v (x,r ) and some r < min£u _1 (disty(y, F(x))), then 
||x — x||y > r for all x G F~ 1 (y) fl B v {x, R). 

A variant of part a) for finite-dimensional Hilbert spaces has been proved 
in 0. This variant will be used in the proof of Theorem [9] together with 
a Galerkin approximation. For this we let (wk)kLi denote an orthonor¬ 
mal basis of V, and for N E N we consider the finite-dimensional subspace 
Vn spanjty!, ... ,Wn} C V. The orthogonal projection from V onto Vn 
is denoted P^. We will then use the following reformulation of pj. Theorem 

3-11- 

Theorem 10. Letx,y G Vn, R > 0 and l < 0, let Fn '■ By N (x,R) —> CC{Vn) 
be a multivalued mapping, and lety G F N (x). If F N is use, if\\y — y\\v < ~IR 
and if for every x G B Vn (x, R) there exists some y G F N (x) such that 

(y-y,x- x) v < l\\x - x\\ 2 v , 

then there exists some x G B Vn (x c , — ^j\\y — y\\v) with x c := x + i^(y— y) and 

y e f n (x). 
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The remainder of this section is devoted to the 
Proof of Theorem [3 

Statement a), special case: We begin with the case x — 0, y — 0, in which 
\\y\\v < —IR, and dehne yjv '■= PnV■ Clearly, B Vn (0,R ) C L>y(0,i?) C 
V and 11 yjy\\v < —IR, and the mapping F N : B Vn (0,R ) =1 Vn given by 
Fn(x) PnF(x) is well-defined with y/v G Fjv(O). Moreover, Fn satisfies 
the assumptions of Theorem flOl on By N ( 0, R) with data xn = 0, y^ = 0 and 
Vn- 

i) The mapping Fjy is bounded on By N (0, R) with convex and compact 
values: As Vn is finite-dimensional and Pn G £(V,Vn), this follows 
from Lemma [TJ 

ii) The mapping Fn satisfies an ROSL-type condition: Let x G By N (0, R). 
By assumption, there exists some y G F(x) such that (y—y, x) < /||a;||y, 
which implies that 

(Rnv - vn, x)v = (y-y, Pnx)v = (y-y, x) v < l\\x\\y. 

iii) The mapping Fn is use: Assume that Fn is not use at x G By N (Q,R). 
Then there exist s > 0 and two sequences (xk)ke¥i C By N (0,R) and 
(Vk)ke¥i C V N with lim^oo x k = x and y k G F N (x k ) such that 

disty(yfc, Fn(x)) > e. 

As a consequence of m Theorem 13.1] there exist (vk)ke n C Vn such 
that ||ufc||v = 1 and 

(■ Vk,yk)v > <Tv N (v k , B Vn (F N (x), e)) = <Jy N (v k , F N (x)) + £. 

As dim V N < oo, there exists v G Vn such that \\vk — v\\y —> 0 along a 
subsequence. Hence 

(v, yk)v — (v - v k , y k )v + (vk, Vk)v 

>- sup \\F N (x')\\v\\v - V k \\v + ay N (v k ,F N (x)) + £ / 4 N 

x'cb Vn (o,r) 

<Tv n (v,F N (x )) +£ 
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as k —* oo, because sup x , 6Bv , ||F A r(x , )||y < oo according to a). On 
the other hand, since F is c-uhc, F is uhc, and for any v G Vat, we have 

limsup <Jv N {v : F N {xk)) — limsup sup (PnV, v)v — limsup sup (y,v)v 

k—>oo k —^oo y£F(xk) /c—>-oo y£F(xk) 

= lim sup ay(n, F(xfc)) < a v {v,F{x)) 

k—yoo 

= sup (y,v)v= sup (P N y,v)v = ct Vn (v,F n (x)). 

y&F(x) y&F(x) 

This contradicts (J4J , and hence Fn is uhc. 

Therefore, Theorem flOl yields for every N G N some x n G B Vn (—yiVn, ~ Jfll^Jvllv) 
with 0 G F n (x n ). Rewrite x N = -^(vn + vn) with ||vjv||v < WvnWv < \\y\\v- 
Then there exists some v G V with ||u||y < \\y\\v and such that Vjy —*■ v along 
a subsequence N' C N. Since moreover y^ —>■ y as N —>■ oo, we infer that 

x N x := -^(y + v) G B v (-±y, -^||y||y) as N' 3 N ->■ oo. 

Furthermore, since 0 G Fn(xn), there exist elements </?tv G F(xn) with 
Pn^n = 0 for iV G N, which implies that (<Pn, Wk) —>• 0 as N —> oo for every 
k G N. Since (tffc)fc is an orthonormal basis of X and the sequence, (^n)n is 
bounded as F is bounded, it follows that tp N —^ 0 as N —)■ oo. For arbitrary 
v G V", we thus find that 

0 = lim < lim sup oy (v, F(x N )) < o v (v,F(x)), 

N^oo 7V->-oo 

because F is c-uhc. This implies 0 G F(x). 

Statement a), general case: Consider x G V, y G V and the map G : 

By{ 0, R) —» CBCiy ) given by 

G(z) := F(z + x) -y, 

Clearly y 0 y — y G F{x) — y = G(0). For any x G B v (x, R ), there exists 
0 G By({), R) such that z + x = x, and by assumption, there exists some 
y G F(x) such that 

(y-y,x- x) v < l\\x - x\\ 2 v = l\\z\\ 2 v . 

But then y' := y — y G F{x) — y — G(z) satisfies 


(y'-yo,z)v = (( y-y ) - (y-y),x-x) v = (y -y,x -x) v < l\\z\\ 2 v , 


so that G satisfies all assumptions of Step 1, which guarantees the existence 
of some z 0 £ By( 0, R ) with 0 £ G(z 0 ) and 

Iko + YiV oik < — ^7112/o 11 v • 

Setting x := x + zq we obtain y £ F{x) and 

Ik - x c\W = \\x-x- ±(y - y)\\ v = |ko + jiVolW < -^IMk = ~i\\y ~ v Ik- 
Statement b): Assume that y £ F(x) for some x £ By(x,r). Then 
dist y(y, F(x)) < ey(F(x), F(x)) < co{\\x — x||v) 


implies 


x — x\\y > milieu 1 (disty(y, F(x))) > r. 


□ 


4 A reformulation for Gelfand triples 

The aim of this section is to adapt the above solvability theorem to a situ¬ 
ation in which the multivalued operator consists of two parts with different 
properties. Theorem [TT] improves the approach presented in [4] by consider¬ 
ing the problem in a space that is adapted to the composed operator. We 
postpone a comparison of both results to Remark [13] at the end of this sec¬ 
tion. The most prominent setting, in which such a splitting occurs, will be 
discussed in the extended example in Section [7] 

Let (V, || • ||y, (•, -)y) and (H, || • j \h, (■, • )h ) be separable Hilbert spaces 
such that V is densely and continuously embedded into Ft with embedding 
constant Cyn > 0. Identifying H with its dual H*, we then have embeddings 


Here i* denotes the dual of i, and this map is injective due to the density of 
i{V) in H. As usual, we regard V as a subspace of H and H as a subspace 
of V*, writing simply v £ Ft instead of i{y) and w £ V* instead of i*(w) for 
v £ V , w £ H. With these simplifications, we have 

(y, x) = (y, x) H Vieh, y £ H. 
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In the following, we fix constants 

ly < 0 and Ih < ~W/ c vh■ (5) 

Then the bilinear form 

(xi,x 2 ) (xi, x 2 )w '■= -ly(x i,x 2 ) v - Ih(xi,x 2 )h, x 1) x 2 eV 

is a scalar product which induces an equivalent norm || • \\w on V. In the 
following, B w (x , R) denotes the ball with radius R > 0 w.r.t. || • \\w centered 
at x G V. Moreover, for y G V*, we denote by ||y||w* the dual norm induced 

by II • || w, he., ||y||w* = sup f° r V *= F*. We also denote by Jw '■ 

xgv\(o} l|3:|lw 

V* —> V the corresponding canonical isometric isomorphism given by 

(J w <p,v)w = {<p,v) Vu G V, (f G V*. 

The following theorem is a variant of Theorem [9] for composite operators. 

Theorem 11. Suppose that lyjn £ K. satisfy (EP, and, let x G V , y G V* 
and R > 0. Moreover, let 


F v : B w (x, R) ChG CBCiV*) and F H : B w (x, R) ChG CBC(H) 


be bounded and c-uhc, and let F : B\y(x, R) =4 V* be given by F = Fy + Fh, 

i.e. 

F{v ) := {y v + y H : Vv e F v (v),y H G -Fff(u)} for v G V. 

Suppose furthermore that there exists yy G F v (x), yn G F H (x) such that 


\\y-y\\w* < R 

with y yv + Vh and 

Mx G B w (x,R) ( (y v -y v ,x - x) <l v \\x - x\\y] 

3 y v G F v (x),y H G F H (x) J | (2/h ~y H ,x~ x) H < l H \\x - x\\ 2 H . 

Finally, let x c := x — \ ~ v)■ Then there exists some 

x G Fh/(x c , |||y - y||w*) satisfying y G F(x). (7) 
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Remark 12. (i) In the case where the embedding of V in H is compact, it 
suffices to assume that F H : B w {x , R) C H —$■ CBC(H) is bounded and uhc, 
because then Fh ■ Bw(x, R) C V —> CBC(H) is bounded and c-uhc. 

(ii) The assumption (J6]) arises naturally in applications, and it is the reason 
for using the mixed norm || ■ \\w which then gives rise to optimal estimates, 
as explained in (iii) below. Nevertheless, sufficient assumptions can easily be 
formulated in terms of || ■ ||y and || • ||# by using the estimates 



for x G V and 



for y E 1/*. 


(iii) In the case where Ih G [0, — -ty-), (0) and (jHJ) imply the estimate 

C VH 




— ly ||a: — x , 


— Ir\\x — x t 


and therefore 



In the case where Ih <0, (0) and (JSD imply the estimate 



(9) 




( 10 ) 


Hence a negative one-sided Lipschitz constant Ih of Fh improves the estimate 
for ||x — x c ||zf, which is important in the case where Ih W < 0. 

Proof of Theorem{Tll We apply Theorem [9] to the Hilbert space (V, (■, -)w), 
the map F 0 := Jw ° F : Bw{x, R) =4 V in place of F and yo := JwV , Vo '■= 
JwV £ V in place of y , y, respectively. We check that F 0 satisfies the as¬ 
sumptions of Theorem 0 
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i) One-sided Lipschitz property. Let x G B w (x,R). By (E]), there exist 
Dv G Fy(x) and yu G F H (x) such that, denoting y := JwilJv + Vh), we 
obtain y G F 0 (x) and 

(■y - yo,x- x) w = (y v + y H -y,x- x) 

= (yv-yv, X -x) + (y H -VH, X- x)h (n) 
< ly\\x — x\\ v + l H\\X — X\\ H = —\\x — x\\ w 

Therefore, F is relaxed one-sided Lipschitz relative to x and y on 
B w (x,R) with constant / = —1 w.r.t || • \\w- 

ii) Properties of the images. For any x G B w (x,R), we have Fy{x) G 
CBC{V*) and F H (x) G CBC(H). Consequently, F H (x) G CBC(V*) by 
Lemma [3 Then Lemma [8] guarantees that F(x) G CBC(V*). Again by 
Lemma 3 it follows that F 0 (x) = Jw{F(x)) G CBC(V). 

iii) F 0 is bounded and c-uhc. The boundedness of F 0 is an easy consequence 
of the boundedness of the maps Fy and Fh- To show that F 0 is c-uhc, 
let {x n ) n C B w (x, R) satisfy x n —x G B w (x, R), and let uGh Using 
that Fh is c-uhc as a map to H by assumption, we then find that 

lim sup cfy(v, Fnixrf)) — lim sup sup (y,v) 

n->oo n^-oo y&F H (x n ) 

= lim sup sup (y, v)h = lim sup a H (v, F H (x n )) 

n—Hxi y£F H (x n ) rwoo 

< a H (v,F H (x)) = sup (y,v) H 

y&F H (x) 

= sup (y,v) = (Ty(v,F H (x)). 

y&F H (x) 

Combining this with the fact that Fy : B w (x,R ) C V —> CBC(V*) is 
c-uhc by assumption, we find that 

lim sup ay(v, F 0 (x n )) = lim sup cry(v, F(x n )) 

n—>oo n—>oo 

= lim sup <Jy(v, F v (x n ) + F H (x n )) 

n —^oo 

= limsup{cry(u, F v (x n )) + af(v, F H (x n ))} 

n—too 

< lim sup cry(v, F v (x n )) + lim sup o~y(v, F H (x n )) 

n—> oo n—> oo 

< <7y(v, Fy(x)) + <Ty(v, F H (x)) = <7y(v, Fy(x) + F H (x)) 

= <Xy(v, F(x)) = (Ty(v, F 0 (x)), 
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and thus F 0 is c-uhc. Note that oy is defined here w.r.t. (•, -)w- 

As a consequence, Theorem [9] applies with / = — 1 and yields the desired 
statement. □ 

Remark 13. In . elliptic partial differential inclusions with ROSL right- 
hand sides have been considered as ROSL operator inclusions. Error esti¬ 
mates for Galerkin approximations have been obtained directly without usage 
of the mixed norm || • \\ w . In that case, the core of such an estimate is an 
inequality similar to CD, but of the shape 

(y - y,x - x) v < l v \\x - x\\y + l H \\x - x\\ 2 H 

< ly\\x — x\\y + max{0, Ih}c\,' H ||x — x\\y, 

because the fixed point argument must be carried out in V where the differ¬ 
ential operator is defined. Therefore, the ROSL constant of Fh cannot be 
exploited when Ih < 0, i.e. when this is most desirable. In this situation, 
considering the inclusion in V equipped with the mixed norm || • ||u/ as above 
yields estimates @2) and (ITU . 

5 Inverses of ROSL mappings 

The properties of relaxed one-sided Lipschitz mappings have been studied, 
e.g., in [lO], [nj and other works of the same author. At that time, no 
quantitative information about solutions of algebraic inclusions was available, 
and thus only qualitaive properties of these mappings and their inverses could 
be given. With Theorem [9] at our disposal, we can now prove some basic 
properties of the inverses. 

Theorem 14. Let F : V —> CBC(V) be c-uhc, bounded on bounded sets and 
l-relaxed one-sided Lipschitz with l < 0. Then its inverse F _1 : V =1 V has 
nonempty weakly sequentially closed images, it is —j-Lipschitz, and it is 0- 
relaxed one-sided Lipschitz. Moreover, the images of F^ 1 satisfy the explicit 
and implicit bounds 

\\F-\y)\\v<-nm\\v-j\\y\\v, ( 12 ) 

diamy (F" 1 (y)) < —y sup dianiy F{x) ( 13 ) 

xGF~ 1 (y) 

for any y G V. In particular, F" 1 is bounded on bounded sets. 
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Proof. Properties of the images of F 1 . For arbitrary y G V, apply Theorem 
E to the data F, x = 0, y = y and arbitrary y G F(0) to find F 1 (y) ^ 0. 

Let y G V, x G F~ 1 (y) and x' G V. By the relaxed one-sided Lipschitz 
property, there exists some y' G F(x') such that 

-|| yf ~ y\\v\W - x\\ v < (y' -y,x' - x)y < Z||a/ - x||^ 

and thus 

\\x'-x\\ v <-)\\y'-y\\ v . (14) 

Considering x' = 0, we conclude that 

IMk <-W -y\W< -ih'Wv - mv < -}F(o)lk - ylblk, 

so that bound (TT21) holds. Moreover, considering x' G F” 1 (?/), we deduce 
from inequality (ITT)) that 

I w - x \\v < - jll y' - y\\v < -J dianiy F(x'), 

so that estimate (TT3|) holds. 

Let y G V and x G V, and let (xk)k&ti C F" 1 (|/) be any sequence such 
that Xk — x as k — * oo. As y G F(xk), we have (y,v) < ay(v, F(x/-)) for all 
v G V and k G N, and since F is c-uhc, it follows that 

(y,v) v < lim sup a v (v, F(x k )) < a v (v,F(x)) 

k —^oo 

for all v G V, which shows y G F(x) and hence x G F" 1 (?/). Therefore, 
F~ l (y) is weakly sequentially closed. 

One-sided and Lipschitz estimates. Let y,y' G V and x G F" 1 (?/). By 
Theorem El there exists some x' G F~ 1 (y / ) such that 

x' G B v [x + i(j/' - y),- yi\W ~ y\\v)- 

In particular, 

\W -AW < - jh' -y\W, 

so that F” 1 is —j-Lipschitz. Writing x' = x + h{y' — y) + v with v G 
B v{0 ,-yM - y\\v), we find 

W - V, X* - x)v = {y' - y, ±(y’ - y) + v) v 

< MvW - y\w + hh' ~ yfv - °> 

so that F~ l is 0-relaxed one-sided Lipschitz. □ 
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Without additional structure, it seems difficult to say more about the 
properties of the inverse. Some multivalued mappings arising in control 
theory or from uncertainties are explicitly given in parameterized form and 
allow a more detailed analysis. The existence of such a parametrization 
implies weak-strong continuity of the multifunction and hence is a substan¬ 
tially stronger assumption than the compact upper hemicontinuity required 
in Theorem El 

Proposition 15. Let (U,du) be a metric space, and let F : V —$■ CBCiV ) be 
parameterized by a function f : V x U —> V satisfying 

a) F(x) = U u£ uf{x,u) for all x E V, 

b) w h->- f(x,u ) is continuous for all x E V, and 

c) x i—^ f{x, u ) is continuous from V endowed with the weak topology to V 
endowed with the norm topology, bounded on bounded sets, and l-one¬ 
sided Lipschitz for all u E U in the sense that 

(f(x, u ) — f{x ', u), x — x')v < l\\x — x'\\y for all u E U 

with a constant l < 0 that is independent of u. 

Then F” 1 : V =4 V is parameterized by a function g : V x U —> V satisfying 

d) F~ l iy) = U ueu g(y,u) for all y EV, 

e) u i —y giy, u ) is continuous for all y E V, and 

f) y I—>■ g(y,u ) is —y -Lipschitz and 0-one-sided Lipschitz for all u EU . 

In particular, if U is connected or path connected, the images of F and F _1 
inherit these properties. 

Proof. Applying Theorem fill to the functions x f(x,u ), u EU, yields the 

existence of inverses g(-,u ) : V =$ V, u EU, given by 

g(y,u) ■= {x E V : f(x,u) = y}, 

that are well-defined, —y-Lipschitz and O-ROSL, so that f) holds. If y E V, 
u EU and x,x' E g(y,u), then 

0 = if (x, u) - f{x' , u), x' - x) v < l\\x - x'\\y, 


15 


which enforces x = x', so that g is a single-valued function. 

Let y G V and x G F~ l (y). Then y G F(x), and hence there exists u G L 
with y = f(x,u), so that x = g(f(x,u),u ) = g(y,u), so that d) is valid. 

Let y G V and u G U be arbitrary, and let (u n ) n C U be such that 
du(u,u n ) —> 0 as n —> oo. By Lipschitz continuity of y (->■ g(y,u n ), we have 

\\a(y,u) - g(y,Un)\\v < —yII f(g(y,u),u n ) - f(g(y,u n ),u n )\\ v 
< —jII f(g(y,u),un) - f(g(y,u),u)\\ v - j|| f(g(y,u),u) - f(g(y,u n ),u n )\\ v 
= —jll f(g(y,u),u n ) - f(g(y,u),u)\\v -> 0 as n -> 00 . 

Hence u g(y,u) is continuous, and we have verified e). □ 

6 Numerical solution of algebraic inclusions 

We propose the algorithm given below in (Il5p for the computation of a so¬ 
lution of the algebraic inclusion y G F(x), where y G V is given and F is 
Z-ROSL and L-Lipschitz with l < 0 and a moderate Lipschitz constant L > 0. 
According to estimate (fI7j) below considered for n — 0, the method, when 
analyzed without round-off errors, i.e. with £„ = 0 for all n G N, finds a 
solution x with 

11% ~ X\\ V < 2 Z+L dist ^(f/, F(x 0 )). 

When applied in this context to the points x = x 0 and y = Proj^(y, F(x 0 )), 
Theorem |9] guarantees the existence of a solution x! of the same inclusion 
with 

x' G B v (x c , -jjdist v (y,F(x 0 ))) and x c = x 0 + ±(y - Proj v (?/, F(x 0 ))), 
which means that 


Iko - x'\\v < -y dist v (y, F(x 0 )). 

Therefore, the solution x of the numerical method is up to a factor 2 +lJJ as 
close to the initial value Xo as the theoretical estimate. This is not necessarily 
true for an arbitrary scheme in the set-valued context, and therefore an 
interesting feature of this algorithm. 

Moreover, it is currently unclear whether a continuous and Z-ROSL multi¬ 
valued mapping possesses a selection or a parametrization that is continuous 
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and (uniformly) /-one-sided Lipschitz. This means that, in general, it is im¬ 
possible to apply a standard numerical method to a single-valued selection 
/ of F and to compute in this way a solution of the multivalued problem. 
The most promising construction in this direction has been published in [2], 
where set-valued mappings were parameterized by generalized Steiner points 
of their images. 

The basic technique behind the following proposition is the same as in [51 
Proposition 4.1]. There are, however, some additional difficulties, because 
the images of F” 1 are not compact and iterates cannot be computed exactly 
in the current setting. Computational errors will therefore be modelled by a 
sequence (£ n ) ra C V. 

Proposition 16. Let y E V, and let F : V —>■ CBC{V) be c-uhc, L-Lipschitz 
and l-relaxed one-sided Lipschitz with l < 0 such that 0</c:=—4<1. For 
arbitrary xq E V, define (v n ) n C V and (x n ) n C V by 

v n := V - Proj v (y, F(x n )), x n+ i := x n + F Vn + f n , n E N, (15) 

where (£ n ) n C V is an arbitrary sequence such that ||£n||v < oo. Then 

the sequence (y n ) n C R + given by 

n 

Vn • 'y ^ ^ HCn—fc||v 
k =0 

satisfies Yl^LoVn < oo, and the sequence (x n ) n converges to some x E F” 1 (?/) 
with estimates 

disty(x ri , F~ l (y)) < -^Huollv + Vn- i, (16) 

OO 

\\x n -x\\v < -%£;\\vo\\v + '52vj- ( 17 ) 

j—n 

Proof. According to Lemma H the projection Projy(r/, F(x)) is a singleton 
for every x E V, so that the sequences (v n ) n and (x n ) n are well-defined. Note 
that 

OO OO OO 

J^vj = < 00 

j=o j =o j =o 

by the Cauchy product formula. Applying Theorem [U] for every n E N with 
x = x n and y = Proj v (y, F(x n )), we find x n E F~ l (y) such that 

disty(x n+ i, F^ 1 (y)) < \\x n+1 - x n \\ v 

||%n T 2 l^ n F "^nllv A 2/|l^n||v T ||Cn||v- 
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By TheoremdH the preimage F 1 (y) is weakly sequentially closed, and there¬ 
fore, by Lemma El there exist points x n G F _1 (|/) such that 

\\x n+1 -x n \\v = disty(x n+ i,F _1 (|/)) Vn G N. 

It follows from inequality (fT8li that 

||w n+ i||y = disty(r/,F(x n+ i)) < e v (F(x n ),F(x n+ i)) < L\\x n - x n+ i\\ v 

= L disty(x n+ i, F“ 1 (|/)) < -7§||Un||y + £||£n||y = K IM|v + L\\€n\\v, 

so that 

||u n ||v < K n |ko||y + Lrjn- 1 , 
and, again because of we have 

disty(x n , F~ l (y)) < ~Yi\\ v n-i\\v + ||£n-i||v 

< — !l W~ il^ollv + I^Vn-2 + ||£n-l||v = ~ \ \ Vq \ \ V + T/ n -l, 

which is fflfijl . Then 

H-^n+l *^n||y — 27 11 11W T ||UV 

< -%\\v 0 \\v + KVn -1 + 11 £,n 11 V < + Vn 

implies that for any n > k, we have 


n —1 


n—1 


\x n ~ X k \\ V < ^ ll^'+l - X j\W < ^(-fr 
j=k j=k 


V 


+ Vj) 


( 19 ) 


< - 


^ -YiU v o\WT^ + “>• 0 as k -7 oo. 

j=k 


In particular, the sequence (x n ) n C V is Cauchy and hence converges to some 
x G V. Since F” 1 (r/) is weakly sequentially closed and inequality (TT6|) holds, 
Lemma El guarantees that x G F~ l (y). Estimate (TT71) follows from (fTUj) by 
passing to the limit n —> oo. □ 


7 Example 

In this section, we consider a class of systems of elliptic differential inclusions. 
Scalar partial differential inclusions with ROSL right-hand sides have been 
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studied in [3]. Existence and relaxation theorems have been proved in a 
more general context. For a recent contribution, we refer to [?]. Elliptic 
partial differential inclusions with multivalued mappings given in terms of 
sub differentials have been studied, e.g., in the monograph |6j. 


7.1 A system of elliptic differential inclusions 

We consider the system 

(—Aui(x),—Au 2 {x)) E f(x,ui(x),u 2 {x)), x E fl, 
u\{x) = u 2 {x) = 0 , 16 


( 20 ) 


of elliptic partial differential inclusions, where C R d is a bounded domain, 
and let / : x R 2 —y CC(R 2 ) be a multivalued mapping with the following 
properties. 

Al) The mapping / is Caratheodory in the sense that x i —> f(x,s\,s 2 ) is 
measurable for any (si,S 2 ) G R 2 and (si,s 2 ) f(x,s\,s 2 ) is continu¬ 

ous for almost every x E fh 

A2) The mapping / is uniformly Z/-ROSL in the sense that for almost every 
x G 12 and every s = (si, s 2 ) G R 2 , t = (ti,t 2 ) G R 2 and 77 = ( 771 , 772 ) G 
f(x,si,s 2 ), there exists p = (pi, p 2 ) G f(x,ti,t 2 ) with 

(. P~r),t-s) < lf\\t — sHi, 

where (•, •) denotes the Euclidean scalar product on R 2 . 

A3) The mapping / is linearly bounded in the sense that there exist a G 
L 2 {yi) and /? > 0 with 

||/(x,si,s 2 )||2 < a(x) +^||(si,s 2 )||2 for a.e. x £ Q, V(si,s 2 )gR 2 . 


The weak formulation of (l20l) is as follows. A pair (ui,u 2 ) G Hq(Q) x Hq(Q) 
of functions is called a weak solution of (l20li if 


3 hi, h 2 G L 2 (Q) s.t. 

(Vui, Vw) L 2 = (hi, w) L 2 \/w G 
(Vu 2 , Vw) L 2 = (h 2 ,w) L 2 \/w G Hq(Q) 

(hi(x), h 2 {x)) G f(x,Ui(x),u 2 (x)) for a.e. x E 


( 21 ) 
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To simplify the notation, we denote V := x Hq(Q) and H := L 2 (fi) x 

L 2 (Q) with V* = if _1 (0) x i/ _1 (h2). The spaces V and H are equipped with 
the scalar products 

(u,v) v ■= («i,ui)hi + (u 2 ,v 2 ) ( h,g) H := + (, h 2 ,g 2 ) L 2 

and the corresponding norms. The duality pairing between 1/ and 1/* is given 
by 

(ip,u) := + (ip 2 ,u 2 ). 

Note that V C H C Id* is a Gelfand triple. It can be shown that the 
set-valued Nemytskii operator given by 

Nf(u) {h G H : h(x) G f(x,u(x)) a.e.} 

is a continuous mapping Nf : H —>■ CBC(Il), which is also //-ROSL. If, 
in addition, the mapping (si,S 2 ) l— t f(x,Si,S 2 ) is L /-Lipschitz w.r.t. the 
Euclidean norm for all x G hi, then Nf is L/-Lipschitz as well. 

We can rewrite (l2l| as an operator inclusion 

0 G (Aiti, Au 2 ) + Nf(ui,u 2 ). (22) 

To comply with the notation in Section 01 we denote F v := (A, A) : V —* V* 
and Fh := Nf : H —^ CBC(H) with one-sided Lipschitz constants ly — — 1 
and Ih = If- As the embedding V C H is compact, the continuity of iVj 
is, according to Remark fl2lfi). sufficient for the application of Theorem [9] 
provided that If < l/cy H . We stress that, by definition of the embedding 
constant Cyn, the quantity 1 /cy H is simply the first Dirichlet eigenvalue of 
—A on fh We have seen in the proof of Theorem fill that the norm 

\\ u \\w = Mv ~ l f\\ u \\H for u G V 

captures the one-sided properties of the composed mapping F v + F H : V —> 
V* in an optimal way in the sense that it is /-ROSL with / := —1 w.r.t. || ■ \\w 
We find that 


e w *(F H (u ), F h {u)) < 

Lf 


< , 

— / 1 


~—Ih 


U — U\\h < —f 


l H 

Lf 

~—h. 


■e H (Nf(u),Nf(u)) 


I — II ijif || — || 

\u-u\\ w = 1 \\u -u\\ w 
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for all u, u G R. In order to check the assumptions for the iterative algorithm 
m , we distinguish two cases. 

Case 1: If < 0. In this case Fy : (R, || • ||w) —¥ (R*, || • ||w*) is 1-Lipschitz, 
and thus Fy A Fh : (R, || • ||w) -A (V*, || • ||w*) is L-Lipschitz with L := 
1 + • Thus we can compute a solution of system (1201) or, equivalently, 

1 C VH l f 

operator inclusion (122|) by applying the iterative algorithm (fT5]h provided 
that L < —21 = 2, or, equivalently, 


Case 2: If G [0,In this case it follows from the estimates in Re- 
mark R^Tii) that Fy : (V, II • II w) —>■ (R*, II • ||w*) is l , -Lipschitz, and 

1 c VH l f 

thus Fy A Fh ■ {V, || • ||w) —> (R*, || ■ ||w*) is 1 ^Y i aLf -Lipschitz. As a conse- 

1 c VH l f 

quence, the iterative algorithm (1T5|) applies in this case if 

Lt < -ji-2 lt : 

1 C VH 

which a posteriori requires If < ^ 3 — since If < Lf. 

7.2 Computational considerations 

We are looking for a solution w G R to the operator inclusion (l22]h Given any 
initial value Uq G R, the numerical routine proposed in (TT5h consecutively 
constructs a sequence (u n ) n C R of approximate solutions that converge to 
a solution of (l22lh In the present context, the iteration reads 

U n -\~\ U n 21 Jw P^oj W* ( 0 ) Fy(u n ) A Fh(Uti)') 

2 1 J\\' Fy [Un) A J\y Nf^Un ')'), (23) 

because Jw ■ (R *, II ' II w*) > (R, II • || w) is an isometrical isomorphism. 
From a computational perspective, it may be advantageous to recast the 
optimization problem 

|| JwFy{u n ) A Jwh\\w = min! subject to h G Nf(u n ) 
with pointwise inequality constraints into an unconstrained dual problem. 
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Lemma 17. Let X be a Hilbert space with inner product (•, -)x and norm 
|| • ||x? and let A G CBC(X). Then the optimization problems 

±||x|ft = min! subject to x G A (24) 

and 

|||x||^ + ctx(— x, A) = min! (25) 

on the entire space X possess the same unique solution. 

Proof. It is well-known (see, e.g., [8, Proposition 7.4]) that (T24li possesses a 
unique solution x* G A, which is also the unique solution of the variational 
inequality 

(x, a — x)x >0 Vagi (26) 

in A. As the function x hg |||x||^- + ax(—x, A) is convex, lower semicontinu- 
ous and coercive, problem (1251) possesses at least one solution. Any solution 
x* G X of (125]) satisfies the necessary optimality condition 

0 G x* - da x (—x *, A), 

where dax{-, •) denotes the sub differential w.r.t. the first variable. This im¬ 
plies 

x* G d<Jx(— x*, A) = argmax aeA (—x*, a)x C A, (27) 

so that x* G A. Moreover, (1271) implies 

(—x*,x*)x > (—x*,a)x for all a G A, 

so that x* solves (1251) and hence x* = x*. □ 

In the situation of our example, the support function in the dual problem 
(125|) can be computed explicitly. 

Lemma 18. For all u, v G V we have 


a w (v,J w N f (u))= / a(v(x), /(x, u(x))dx, 


IQ 


where cr(-, •) denotes the support function on R 2 . 
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Proof. Given u,v G V, we construct some h v G Nf(u ) such that 

(v, h v ) H = max (v, h) H . 


heN f {u) 


( 28 ) 


Theorem 8.2.11 in [1] on marginal maps ensures that the multivalued map¬ 
ping G : Q —» C£>C(R 2 ) given by 

G(x) := {t G f(x, u(x)) : (v(x),t) = max (u(a;), s)} 7 ^ 0 

s€f(x,u(x)) 

is measurable, and [H Theorem 8.1.3] ensures that G possesses a measurable 
selection h v : hi — > R 2 , i.e. 

h v {x) G G(x) C f(x,u(x)) for a.e. x G hi. 

The linear growth bound A3) ensures that h v G H, and therefore h v G 
Nf(u n ). By monotonicity of the integral and the construction of h v , we have 

(v, h v ) H ~ (v, h) H = / (n(x), h v (x) - h(x))dx > 0 


for all h G Nf(u), so that h v satisfies condition (l28j) . By construction of h v , 
we hnd 


a w (v, J w Nf(u)) — sup ( v,J w h) w = sup (v,h) H 

h£Nf(u) h£Nf(u) 



□ 


The following statement is a consequence of Lemmas |TT] and fl8l 

Corollary 19. In the iteration defined by [23) . the function —2 l(u n+ i—u n ) G 
V is the unique minimizer of the functional I : V —> R given by 
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7.3 Some numerical results 


We first consider the problem 

4 Ini 2 

(—Aiti, —A u 2 ) G —— • . — r^u + IfU + B R ( 0) (29) 

9 1 + |n| z 

on fl = (0,1) with u = 0 on <912 and If < 0. Elementary computations 
show that the right-hand side is If -ROSL and Lj-Lipschitz with a constant 
Lf = \ — lf. According to Section 17.11 the composed mapping 


F v + F H : (V, || • ||w) -;> (Y*, II • II w*) 
is /-ROSL and L-Lipschitz with constants / = —1 and 


Lf _ _ 1 , V2-t / 


= 1 + -ATt 1 < 2, 

7 T Z —lf 


so that algorithm 


L — 1 —1 , / 2 I 

l/ c VH~h ’ K "~ l f 

is applicable with a theoretical speed of convergence 

K :=-L/2i = ±(1 + ggg) < 1. 

The results for parameters If — — 1 and R = 10 and initial values 

uoi{x) = |sin(27rx), U 02 = ^ sin(167nr) 

displayed in Figure Q] show that the bound (fTTH) is realistic in this case. The 
residual 


r n = distiy.(0, A u n + N F (u n )) 

is approximately halved in every iteration, while the theoretical bound guar¬ 
antees a reduction by a factor 


K — \ + 


3 

4(7T 2 -|-1) 


« 0.569. 


We now consider the problem 

(—Aiti, -A u 2 ) E (—U!U 2 + 1 - Mi + x, —uiu 2 + 1 - u 2 + x) + B r ( 0) (30) 

on 0 = (0,1) with u = 0 on dfl. The application of algorithm (1T5|) to this 
problem is a-priori not theoretically justified. Examples not included here 
show that the iteration can indeed diverge. For most nonnegative initial 
values of moderate magnitude, however, the iteration converges as depicted 
in F igure [21 where R = 5 and 

, , , , (x-O.l) 2 , , , (x-0.8) 2 

u 0 i(x) = x(l — x)e , u 02 (x) = x(l — x)e 5 xG[0,l]. 

In this case, the decay of the residuals justifies a-posteriori the use of the 
method and the validity of the result. 
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Figure 1: Iterates of the numerical algorithm (1T5]) applied to inclusion (129|) . 
The residual is measured in the norm || • \\w*- 
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Figure 2: Iterates of the numerical algorithm (TT5l) applied to inclusion (l30]h 
The residual is measured in the norm II ■ ||w». 
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